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Abstract 

Amplitude models are applied to studies of resonance structure in D°—KgK“7t+ 
and —)• KgK''"7t” decays using pp collision data corresponding to an integrated 
luminosity of 3.0 fb~^ collected by the LHCb experiment. Relative magnitude and 
phase information is determined, and coherence factors and related observables are 
computed for both the whole phase space and a restricted region of lOOMeV/c^ 
around the K*(892)^ resonance. Two formulations for the Ktt S'-wave are used, both 
of which give a good description of the data. The ratio of branching fractions R(D° —)• 
KgK+7t“)/H(D° —)• KgK“7T+) is measured to be 0.655 ± 0.004 (stat) ± 0.006 (syst) 
over the full phase space and 0.370 ± 0.003 (stat) ± 0.012 (syst) in the restricted 
region. A search for CP violation is performed using the amplitude models and no 
significant effect is found. Predictions from SU(3) flavor symmetry for K*(892)K 
amplitudes of different charges are compared with the amplitude model results. 
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1 Introduction 


A large variety of physics can be accessed by studying the decay^ D° —)■ KgK 7 t+ and 
K 0 K+ 71 . Analysis of the relative amplitudes of intermediate resonances contributing 
to these decays can help in understanding the behavior of the strong interaction at low 
energies. These modes are also of interest for improving knowledge of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, and CP-violation measurements and mixing 
studies in the D°-D° system. Both modes are singly Cabibbo-suppressed (SCS), with 
the KgK“ 7 r+ hnal state favored by approximately xl.7 with respect to its KgK+Tt” 
counterpart [^. The main classes of Feynman diagrams, and the sub-decays to which they 
contribute, are shown in Fig. 

Flavor symmetries are an important phenomenological tool in the study of hadronic 
decays, and the presence of both charged and neutral K* resonances in each —)■ KgK^Tt^ 
mode allows several tests of SU(3) flavor symmetry to be carried out |^|^. The KgK^Tt^ 
hnal states also provide opportunities to study the incompletely understood Ktt S-wave 
systems [^, and to probe several resonances in the KgK^ decay channels that are poorly 
established. 

An important goal of havor physics is to make a precise determination of the CKM 
unitarity-triangle angle 7 = arg(—Information on this parameteij^ can 
be obtained by studying CP-violating observables in the decays B“ —)■ D , where the 
D° and D° are reconstructed in a set of common hnal states [^|^, such as the modes 
D°—)■ KgK^Tr'*' and D°—)■ KgK’'" 7 r“ |^. Optimum statistical power is achieved by studying 
the dependence of the CP asymmetry on where in three-body phase space the D-meson 
decay occurs, provided that the decay amplitude from the intermediate resonances is 
sufficiently well described. Alternatively, an inclusive analysis may be pursued, as in 
Ref. 10 , with a ‘coherence factor’ [Tl] parameterizing the net ehect of these resonances. 


The coherence factor of these decays has been measured by the CLEO collaboration 
using quantum-correlated D° decays at the open-charm threshold 12 , but it may also be 
calculated from knowledge of the contributing resonances. In both cases, therefore, it is 
valuable to be able to model the variation of the magnitude and phase of the D°-decay 
amplitudes across phase space. 

The search for CP violation in the charm system is motivated by the fact that several 
theories of physics beyond the Standard Model (SM) predict enhancements above the 
very small effects expected in the SM 13 -15 . Singly Cabibbo-suppressed decays provide 


a promising laboratory in which to perform this search for direct CP violation because 
of the signihcant role that loop diagrams play in these processes 16 . Multi-body SCS 


decays, such as D°—KgK”?!'^ and D°—)■ KgK’^Tr”, have in addition the attractive feature 
that the interfering resonances may lead to CP violation in local regions of phase space, 
again motivating a good understanding of the resonant substructure. The same modes 
may also be exploited to perform a mixing measurement, or to probe indirect CP 

violation, either through a time-dependent measurement of the evolution of the phase 


'^The inclusion of charge-conjugate processes is implied, except in the definition of CP asymmetries. 
^Another notation, bs = 7, exists in the literature. 
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Figure 1: SCS classes of diagrams contributing to the decays )• KgK^Tt^. The color-favored 
(tree) diagrams (a) contribute to the Kg^ Kg7t=*= and (aQ 2 , p)^—> channels, while the 

color-suppressed exchange diagrams (b) contribute to the (aQ 2 ,p)^—^ KgK^, Kq° 2 
and Kq 2 Yi~ TV channels. Second-order loop (penguin) diagrams (c) contribute to the 

(ag 2 , p)^ —>• and Kg^ 2 ~^ Kg7t^ channels, and, hnally, OZI-suppressed penguin annihilation 

diagrams (d) contribute to all decay channels. 


space of the decays, or the inclusive KgK Tf*" and KgK+Tt hnal states 17 


In this paper time-integrated amplitude models of these decays are constructed and 
used to test SU(3) flavor symmetry predictions, search for local CP violation, and compute 
coherence factors and associated parameters. In addition, a precise measurement is 
performed of the ratio of branching fractions of the two decays. The data sample is 


obtained from pp collisions corresponding to an integrated luminosity of 3.0 fb~^ collected 
by the LHCb detector 18,19 during 2011 and 2012 at center-of-mass energies ^/s = 7TeV 
and 8 TeV, respectively. The sample contains around one hundred times more signal decays 
than were analyzed in a previous amplitude study of the same modes performed by the 
CLEO collaboration 

The paper is organized as follows. In Sect.[^ the detector, data and simulation samples 
are described, and in Sect. |^the signal selection and backgrounds are discussed. The 
analysis formalism, including the definition of the coherence factor, is presented in Sect. 
The method for choosing the composition of the amplitude models, ht results and their 
systematic uncertainties are described in Sect.[^ The ratio of branching fractions, coherence 
factors, SU(3) flavor symmetry tests and CP violation search results are presented in 
Sect. Finally, conclusions are drawn in Sect. 


2 Detector and simulation 

The LHCb detector is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < 77 < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
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surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, of charged particles with a relative uncertainty 
that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of 
a track to a primary pp interaction vertex (PV), the impact parameter, is measured with a 
resolution of (15 + 29 /pt) where px is the component of the momentum transverse to 
the beam, in GeV/c. Different types of charged hadrons are distinguished using information 
from two ring-imaging Gherenkov (RIGH) detectors. Photons, electrons and hadrons are 
identihed by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. 

The trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, in which all charged particles with 
Pt > 500 (300) MeV/c are reconstructed for 2011 (2012) data. At the hardware trigger 
stage, events are required to have a muon with high px or a hadron, photon or electron with 
high transverse energy in the calorimeters. For hadrons, the transverse energy threshold 
is 3.5 GeV. Two software trigger selections are combined for this analysis. The hrst 
reconstructs the decay chain D*(2010)’''—?■ with D°—)■ h+h“X, where h^ represents 

a pion or a kaon and X refers to any number of additional particles. The charged pion 
originating in the D*(2010)^ decay is referred to as ‘slow’ due to the small Q-value of the 
decay. The second selection fully reconstructs the decay D°—)■ KgK^Tt^, without flavor 
tagging. In both cases at least one charged particle in the decay chain is required to have 
a signihcant impact parameter with respect to any PV. 

In the offline selection, trigger signals are associated with reconstructed particles. 
Selection requirements can therefore be made on the trigger selection itself and on whether 
the decision was due to the signal candidate, other particles produced in the pp collision, 
or both. It is required that the hardware hadronic trigger decision is due to the signal 
candidate, or that the hardware trigger decision is due solely to other particles produced 
in the pp collision. 

Decays Kg —)■ t&tC are reconstructed in two different categories: the hrst involves Kg 
mesons that decay early enough for the pions to be reconstructed in the vertex detector; 
the second contains Kg mesons that decay later such that track segments of the pions 
cannot be formed in the vertex detector. These categories are referred to as long and 
downstream, respectively. The long category has better mass, momentum and vertex 
resolution than the downstream category, and in 2011 was the only category available in 
the software trigger. 

In the simulation, pp collisions are generated using Pythia 
conhguration 


21 


22 . Decays of hadronic particles are described by EvtGen 23 


with a specihc LHGb 
in which 


hnal-state radiation is generated using Photos 24 . The interaction of the generated 


particles with the detector, and its response, are implemented using the Geant4 toolkit 
as described in Ref. 


26 


25 


3 








3 Signal selection and backgrounds 

The offline selection used in this analysis reconstructs the decay chain D*(2010)’^ —?■ 
with D°—)■ KgK^TirF, where the charged pion from the D*(2010)’'' decay tags the flavor 
of the neutral D meson. Candidates are required to pass one of the two software trigger 
selections described in Sect.[^ as well as several offline requirements. These use information 
from the RICH detectors to ensure that the charged kaon is well-identified, which reduces 
the background contribution from the decays )■ Kg 7 r'^ 7 r“ 7 t° and D°—)■ Kg 7 r^|j.+v^. In 
addition the Kg decay vertex is required to be well-separated from the D° decay vertex 
in order to suppress the D° —)■ K“ 7 r'*’ 7 r"*' 7 r“ background, where a 7 r’*" 7 r“ combination is 
close to the Kg mass. D° candidates are required to have decay vertices well-separated 
from any PV, and to be consistent with originating from a PV. This selection suppresses 
the semileptonic and D° —)■ K'Tt+Tt+Tt” backgrounds to negligible levels, while a small 
contribution from D° —)■ Kg7f^7r“7r‘’ remains in the Am = m(KgK7r7tsiow) — m(KgK7t) 
distribution. A kinematic fit is applied to the reconstructed D*(2010)+ decay chain 
to enhance the resolution in m(KgK 7 r), Am and the two-body invariant masses m(KgK), 
m(Kg7t) and m(K7r) that are used to probe the resonant structure of these decays. This fit 
constrains the D*(2010)’'' decay vertex to coincide with the closest PV with respect to the 
D*(2010)’'' candidate, fixes the Kg candidate mass to its nominal value, and is required to 
be of good quality. 

Signal yields and estimates of the various background contributions in the signal window 
are determined using maximum likelihood fits to the m(KgK 7 t) and Am distributions. 
The signal window is defined as the region less than 18MeV/c^ (0.8MeV/c^) from the peak 
value of m(KgK 7 t) (Am), corresponding to approximately three standard deviations of 
each signal distribution. The three categories of interest are: signal decays, mistagged 
background where a correctly reconstructed D° meson is combined with a charged pion that 
incorrectly tags the D° flavor, and a combinatorial background category, which also includes 
a small peaking contribution in Am from the decay D° —)■ Kg 7 r''' 7 r~ 7 r°. These fits use 
candidates in the ranges 139 < Am < 153MeV/c^ and 1.805 < m(KgK 7 r) < 1.925 GeV/c^. 
The sidebands of the m(KgK 7 t) distribution are defined as those parts of the fit range where 
m(KgK 7 r) is more than 30 MeV/c^ from the peak value. The Am (m(KgK 7 r)) distribution in 
the signal region of m(KgK 7 r) (Am) is fitted to determine the D*(2010)^ (D°) yield in the 
two-dimensional signal region [^. The D*(2010)’'' (D°) signal shape in the Am (m(KgK 7 T)) 
distribution is modeled using a Johnson Su (Cruijff [^) function. In the m(KgK 7 t) 
distribution the combinatorial background is modeled with an exponential function, while 
in Am a power law function is used, /(Am; m^,p, P, 5) = ( ^"" )p - y-^( ^"" )^, 
with the parameters p, P and b determined by a fit in the m(KgK 7 t) sidebands. The small 
D°—> Kg 7 f'' 7 r“ 7 T' contribution in the Am distribution is described by a Gaussian function, 
and the component corresponding to D° mesons associated with a random slow pion is 
the sum of an exponential function and a linear term. These fits are shown in Fig. The 
results of the fits are used to determine the yields of interest in the two-dimensional signal 
region. These yields are given in Table for both decay modes, together with the fractions 
of backgrounds. 
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Table 1: Signal yields and estimated background rates in the two-dimensional signal region. The 
larger mistag rate in the )• KgK'’'7t“ mode is due to the different branching fractions for the 
two modes. Only statistical uncertainties are quoted. 


Mistag Combinatorial 

Mode Signal yield background [%] background [%] 


D°^K°K-7t+ 113 290 ±130 0.89 ±0.09 3.04 ±0.14 

D°^K°K+7t- 76 380 ±120 1.93 ±0.16 2.18 ±0.15 


A second kinematic fit that also constrains the D° mass to its known value is performed 
and used for all subsequent parts of this analysis. This fit further improves the resolution 
in the two-body invariant mass coordinates and forces all candidates to lie within the 
kinematically allowed region of the Dalitz plot. The Dalitz plots 31 for data in the 
two-dimensional signal region are shown in Fig. Both decays are dominated by a 
K*(892)^ structure. The K*(892)° is also visible as a destructively interfering contribution 
in the D°—)■ KgK“7r+ mode and the low-m^o^j region of the D°—)■ KgK+yr” mode, while a 

clear excess is seen in the high-m^o^j region. Finally, a veto is applied to candidates close 


to the kinematic boundaries; this is detailed in Sect. 4.3 


4 Analysis formalism 

The dynamics of a decay D°—;■ ABC, where D°, A, B and C are all psendoscalar mesons, 
can be completely described by two variables, where the conventional choice is to use a pair 
of squared invariant masses. This paper will use = m^(Kg7t) and = m?{Kn) as 
this choice highlights the dominant resonant structure of the D° —)■ KgK^Tt^ decay modes. 

4.1 Isobar models for D^—K^K^Tt^ 

The signal isobar models decompose the decay chain into D° —)■ (R—;■ (AB)j)C contri- 
bntions, where R is a resonance with spin J equal to 0, 1 or 2. Resonances with spin 
greater than 2 shonld not contribnte significantly to the D° —)■ KgK^Tr^ decays. The 
corresponding 4-momenta are denoted poo? Ta, Pb and pc- The reconstrncted invariant 
mass of the resonance is denoted tuab, and the nominal mass ttir. The matrix element for 
the D°—)■ KgK^TifF decay is given by 

m\c), ( 1 ) 

R 
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m{K^^Kn) [GeV/c^ ] 


Am [MeV/c^ ] 


Figure 2: Mass (left) and Am (right) distributions for the —)■ KgK~7T+ (top) and —)• 
KgK'’'7t“ (bottom) samples with fit results superimposed. The long-dashed (blue) curve represents 
the D*(2010)^ signal, the dash-dotted (green) curve represents the contribution of real mesons 
combined with incorrect and the dotted (red) curve represents the combined combinatorial 
and D^—)• Kg7f*'7t“7t^ background contribution. The vertical solid lines show the signal region 
boundaries, and the vertical dotted lines show the sideband region boundaries. 


where is the complex amplitude for R and the contributions Aiji from each 

intermediate state are given by 




R/ 


( 2 ) 


where Rj°(p, |po|, dno) and Bf{q, go, d^) are the Blatt-Weisskopf centrifugal barrier factors 
for the production and decay, respectively, of the resonance R 32 . The parameter p (g) 
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Figure 3: Dalitz plots of the D°—)• KgK 7t+ (left) and D°—)> KgK+7T (right) candidates in the 
two-dimensional signal region. 


Table 2: Blatt-Weisskopf centrifugal barrier penetration factors, Bj{q,qQ,d) 



J Bj{q, go, d) 


0 1 


1 

2 


l+igod)'^ 

l+igd)"^ 


9+3{god)'^ + {god)* 
9+3{gd)'^ + (gd)^ 


is the momentum of C (A or B) in the R rest frame, and po (go) is the same quantity 
calculated using the nominal resonance mass, m^. The meson radius parameters are set 
to d -[)0 = 5.0 (GeV/c)“^ and = 1.5 (GeV/c)“^ consistent with the literature [^,33 ; the 


systematic uncertainty due to these choices is discussed in Sect. 5.2, Finally, m^c) 

is the spin factor for a resonance with spin J and Tr is the dynamical function describing 
the resonance R. The functional forms for Bj(g, go, d) are given in Table and those 
for Gj(m7B,^Ac) Table for J = 0,1,2. As the form for fli is antisymmetric in the 
indices A and B, it is necessary to dehne the particle ordering convention used in the 
analysis; this is done in Table The dynamical function Tr chosen depends on the 
resonance R in question. A relativistic Breit-Wigner form is used unless otherwise noted 


ddl(^AB) — 


- m^e) - ’ 


( 3 ) 
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Table 3: Angular distribution factors, il,j{pY )0 +Pc,PB —Pa)- These are expressed in terms of 
the tensors ^ 1 


J ^./(Tdo+Tc,Pb-Pa) 

0 1 

1 (Pdo + P^)T^^(pb - PA)/(GeV/c)2 

2 (Pdo + Pc)(Pdo + Pc)^m^«/3(Pa - Pb)(Pa “ Pb)/(G eV/c)^ 


Table 4: Particle ordering conventions used in this analysis. 


Decay 

A 

B 

c 



Tt 

K 

K° 

KTK*±, 

^ K^Tt^ 

K° 

71 

K 

D°—i- 7r’'(p=‘',a^), 

p^,a^^ K°K± 

K 


71 


where the mass-dependent width is 

rR(mAB) = Tr [i?j (g, go, ^r)] 


\ 2c/+l 

2 mR ( q \ 


(4) 


channel function used to describe the ao(980)^ resonance 12,35-38 


mAB \qoJ 

Several alternative forms are used for specialized cases. The Flatte (34] form is a coupled- 


(5) 


Tk = 


(^ M 980 )± - - ^[PkkPkk + ’ 

where the phase space factor is given by 

1 

~2 

m: 


Pab = -4~\/("^kk “ - {tiia - 

^KK 


( 6 ) 


and the coupling constants and grin are taken from Ref. 35 , hxed in the isobar model 
hts and tabulated in Appendix The Gounaris-Sakurai 39 parameterization is used to 
describe the p(1450)=‘' and p(1700)=‘' states (^ 40-4^, 


Tr = 




(m2 + /(m^^,m2,rp) -imprp(mKK)’ 


KK 


where 


d{mp) = 


dml 


log 


rrip + 2go 

2mK 


+ 


rur 


27rgo 


m^mp 


(7) 


( 8 ) 
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and 


m: 


/(m^j^,mj,rp) = Tp^ {ql + qlh'{ml){m 

% 


Kk)]- (9) 


The parameter mx is taken as the mean of m^o and mK±, and h'{rni) = jg calculated 


from 


drri^ 


2n 2g(m2) 
h[m ) = ——— log 


vrm 


m + 2g(m^ 


in the limit that mx = mx± = 


mxo. 


2mx 

Parameters for the p resonances 


( 10 ) 

p(1450)^ and 

p(1700)^ are taken from Ref. 44 and tabulated in Appendix [A} 

This analysis uses two different parameterizations for the Ktt S'-wave contributions, 
dubbed GLASS and LASS, with different motivations. These forms include both Kq( 1430) 
resonance and nonresonant Ktt S'-wave contributions. The LASS parameterization takes 
the form 


Tr = / 


TflKn 


4 




"iK5(1430)/ Q 


sin(55 -L 5F)e 


iiSs+Sp) 


( 11 ) 


where f{x) = A exp {hix -t- h 2 x‘^ + h^x^) is an empirical real production form factor, and 
the phases are defined by 

2aq ^ mRrR(mx7i) 


tan 5^7 = 


2 -|- arq'^ 


tan55 = 




m 


( 12 ) 


Ktt 


''K5(1430) 

The scattering length a, effective range r and Kq( 1430) mass and width are taken from 
measurements at the LASS experiment and are tabulated in Appendi x [A} With 
the choice f{x) = 1 this form has been used in previous analyses e.g. Refs. [47| - [49| , and if 
6f is additionally set to zero the relativistic S-wave Breit-Wigner form is recovered. The 
Watson theorem states that the phase motion, as a function of Ktt invariant mass, is 
the same in elastic scattering and decay processes, in the absence of final-state interactions 
(he. in the isobar model). Studies of Ktt scattering data indicate that the S'-wave remains 
elastic up to the Kp' threshold (^. The magnitude behavior is not constrained by the 
Watson theorem, which motivates the inclusion of the form factor f(x), but the LASS 
parameterization preserves the phase behavior measured in Ktt scattering. The real form 
factor parameters are allowed to take different values for the neutral and charged Kq( 1430) 
resonances, as the production processes are not the same, but the parameters taken from 
LASS measurements, which specify the phase behavior, are shared between both Ktt 
channels. A transformed set b' = U~^b of the parameters b = ( 61 , 62 ,^ 3 ) are also defined 
for use in the isobar model fit, which is described in detail in Sect. The constant matrix 
U is chosen to minimize fit correlations, and the form factor is normalized to unity at the 
center of the accessible kinematic range, e.g. |(mxg + -|- mRO — ruxi) for the charged 

Ktt S'-wave. 

The GLASS (Generalized LASS) parameterization has been used by several recent 
amplitude analyses, e.g. Refs. [1^37 38], 

Tr = [Fsinidp + (13) 
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where 6f and 6s are defined as before, and -F, and are free parameters in the fit. 
It shonld be noted that this fnnctional form can result in phase behavior significantly 
different to that measured in LASS scattering data when its parameters are allowed to 


vary freely. This is illustrated in Fig. 10 in Sect. 5.3 


4.2 Coherence factor and C!P-even fraction 

The coherence factor and mean strong-phase difference for the multi-body decays 
D —>■ f and D —f quantify the similarity of the two decay structures 
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In the limit 

—?• 1 the matrix elements for the two decays are identical. For D° —)• KgK^Ti^ the 


coherence factor and mean strong-phase difference are defined by lOl 12 




I MKOK+n- ^ 


2 

Kn 


)dm^^n_dm 




2 

Kn 


^KgK+7t-^KgK- 


where 


■^K9K±7t=F — / l-^KgK±7t=F ^Kn) I '^^KSTr'^^Kn) 


(14) 


(15) 


and the integrals are over the entire available phase space. The restricted phase space 
coherence factor is defined analogously but with all integrals restricted to 

an area of phase space close to the K*(892)^ resonance. The restricted area is defined 
by Ref. [T^ as the region where the KgTt^ invariant mass is within lOOMeV/c^ of the 
K*(892)=*" mass. The four observables i?KgK 7 t) ^KgKn) Rk*k and 5 k*k were measured using 
quantum-correlated i|)(3770) —?■ decays by the CLEO collaboration 12 , and the 


coherence was found to be large for both the full and the restricted regions. This analysis 
is not sensitive to the overall phase difference between D°—)■ KgK’^Ti” and D°—)■ KgK’^Ti'. 
However, since it cancels in (5 kok7i ~ <5k*k 5 this combination, as well as RKgKn and i?K*K, 
can be calculated from isobar models and compared to the respective CLEO results. 


An associated parameter that it is interesting to consider is the CP-even fraction 51 , 




KD+|K°K±7r^)P 


|(D+|K0K±7rT)|2 + |(D_|K0K±7tT)|2 
1 
2 


1 -I- 2PK:gK7t COs(5K0K7t) \j ^KgK7r(l + ^KgKn) 


-1 


(16) 


where the CP eigenstates |D±) are given by ^ [|D°) ± |D°)] and i^KgRTi is the ratio of 
branching fractions of the two D° —)■ KgK^Tt^ modes. As stated above, the relative strong 
phase ^KgKTt is not predicted by the amplitude models and requires external input. 


4.3 Efficiency modeling 

The trigger strategy described in Sect. and to a lesser extent the offline selection, 
includes requirements on variables such as the impact parameter and pr of the various 
charged particles correlated with the 2-body invariant masses and There is. 
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Figure 4: Efficiency function used in the isobar model fits, corresponding to the average efficiency 
over the full dataset. The coordinates and are used to highlight the approximate 

symmetry of the efficiency function. The z units are arbitrary. 


therefore, a significant variation in reconstruction efficiency as a function of and 

This efficiency variation is modeled using simulated events generated with a uniform 
distribution in these variables and propagated through the full LHCb detector simulation, 
trigger emulation and offline selection. Weights are applied to the simulated events to 
ensure that various subsamples are present in the correct proportions. These weights 
correct for known discrepancies between the simulation and real data in the relative 
reconstruction efficiency for long and downstream tracks, and take into account the ratios 
of ^/s = 7TeV to a/s = 8TeV and D°—)■ K3K“7t+ to D°—)■ KgK+7r“ simulated events to 
improve the description of the data. The efficiencies of offline selection requirements based 
on information from the RICH detectors are calculated using a data-driven method based 

decays, where D° —)■ K'Tf*'. These 
A non-parametric kernel estimator 


on calibration samples 52 of D*(2010) 




53 


efficiencies are included as additional weights, 
is used to produce a smooth function describing the efficiency variation in 

the isobar model hts. The average model corresponding to the full dataset recorded in 
2011 and 2012, which is used unless otherwise noted, is shown in Fig. Candidates very 
near to the boundary of the allowed kinematic region of the Dalitz plot are excluded, as 
the kinematics in this region lead to variations in efficiency that are difficult to model. It 
is required that max(| cos( 6 *k 07 i) |, | cos(6'n;K)|, I cos(6'KKg)|) < 0.98, where ^ab is the angle 
between the A and B momenta in the AC rest frame. This criterion removes 5% of the 
candidates. The simulated events are also used to verify that the resolution in 

and is around 0.004 GeV/c^, corresponding to Ci(2MeV/c^) resolution in m(K37r^;. 
Although this is not explicitly accounted for in the isobar model hts, it has a small effect 
which is measurable only on the parameters of the K*(892)^ resonance and is accounted 
for in the systematic uncertainties. 
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4.4 Fit components 


There are three event categories described in Sect. that mnst be treated separately 
in the isobar model hts. The signal and mistagged components are described by terms 
proportional to '^K 7 t)l-^KgK± 7 tT(^K 07 t’''^K 7 t)P) while the combinatorial component 

is described by a smooth fnnction, CK 0 K± 7 iT(’^K° 7 r’obtained by applying to data 
in the m(KgK7r) sidebands the same non-parametric kernel estimator used to model the 
efficiency variation. The same combinatorial background model is used for both D° flavors, 
and the same efficiency function is used for both modes and D° flavors. The overall 
function used in the fit to D°—)■ KgK^TrfF decays is therefore 

-pKgK±7t=F("^Kg7t’= (1 “ /m “ /c)£(’^Kg70 l-^KgK±7tT ("^Kgn’T 

+ m|„) |>lKgK±nT ("iKg„, r (17) 

+ /cCKgK±7rT(”^Kg7t’ 


where the mistagged contribution consists of D°—)■ KgK=*"7r^ decays and /m (/c) denotes 
the mistagged (combinatorial) fraction tabulated in Table 

All parameters except the complex amplitudes are shared between the PDFs 

for both modes and both D° flavors. For the other parameters, Gaussian constraints are 
included unless stated otherwise. The nominal values used in the constraints are tabulated 
in Appendix]^ No constraints are applied for the Ktt F-wave parameters 6i.,3, F, (ps and 
(pF, as these have no suitable nominal values. The Ktt F-wave parameters a and r are 
treated differently in the GLASS and LASS models. In the LASS case these parameters are 
shared between the neutral and charged Ktt channels and a Gaussian constraint to the 
LASS measurements 45 is included. In the GLASS case these are allowed to vary freely 


and take different values for the two channels. 


5 Isobar model fits 


This section summarizes the procedure by which the amplitude models are constructed, 
describes the various systematic uncertainties considered for the models and Anally discusses 
the models and the coherence information that can be calculated from them. 

Amplitude models are fitted using the isobar formalism and an unbinned maximum- 


likelihood method, using the GooFit 54 package to exploit massively-parallel Graphics 


Processing Unit (GPU) architectures. Where y^/bin values are quoted these are simply 
to indicate the fit quality. Statistical uncertainties on derived quantities, such as the 
resonance fit fractions, are calculated using a pseudoexperiment method based on the fit 
covariance matrix. 


5.1 Model composition 

Initially, 15 resonances are considered for inclusion in each of the isobar models: 
K*(892,1410,1680)°’±, K*_2(1430)°’±, ao(980,1450)^, a2(1320)± and p(1450,1700)^. Pre- 
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liminary studies showed that models containing the K*(1680) resonances tend to include 
large interference terms, which are cancelled by other large components. Such hne-tuned 
interference effects are in general unphysical, and are therefore disfavored in the model 
bnilding 36,55 . The K*(1680) resonances are not considered fnrther, and additionally 
the absolnte valne of the snm of interference fractions 56 is reqnired to be less than 


30% in all models. In the absence of the K*(1680) resonances, large interference terms 
are typically generated by the Ktt S'-wave contributions. The requirement on the sum of 
interference fractions, while arbitrary, allows an iterative procedure to be used to search 
for the best amplitude models. This procedure explores a large number of possible starting 
conhgurations and sets of resonances; it begins with the most general models containing 
all 13 resonances and considers progressively simpler conhgurations, trying a large number 
of initial ht conhgurations for each set of resonances, nntil no further improvement in ht 
qnality is fonnd among models simple enongh to satisfy the interference fraction limit. 
Higher valnes of this limit lead to a large number of candidate models with similar ht 
qnality. 

A second procednre iteratively removes resonances from the models if they do not 
signihcantly improve the ht qnality. In this step a resonance mnst improve the valne of 
—21og£, where £ is the likelihood of the full dataset, by at least 16 units in order to 
be retained. Up to this point, the K*(892) mass and width parameters and Ktt S'-wave 
parameters have been allowed to vary in the ht, but mass and width parameters for other 
resonances have been hxed. To improve the qnality of ht fnrther, in a third step, S and 
P-wave resonance parameters are allowed to vary. The tensor resonance parameters are 
known precisely |^, so remain hxed. At this stage, resonances that no longer signihcantly 
improve the ht qnality are removed, with the threshold tightened so that each resonance 
must increase —2 log C by 25 nnits in order to be retained. 

Finally, parameters that are consistent with their nominal values to within Icr are hxed 
to the nominal value. The nominal values used are tabulated in Appendix The entire 
procednre is performed in parallel nsing the GLASS and LASS parameterizations of the 
Ktt S'-wave. The data are fonnd to prefer a solntion where the GLASS parameterization 
of the charged Ktt S'-wave has a poorly constrained degree of freedom. The hnal change 
to the GLASS models is, therefore, to hx the charged Ktt S-wave F parameter in order to 
stabilize the uncertainty calculation for the two corresponding or parameters by reducing 
the correlations among the free parameters. 


5.2 Systematic uncertainties 

Several sources of systematic uncertainty are considered. Those due to experimental issues 
are described hrst, followed by uncertainties related to the amplitnde model formalism. 
Unless otherwise stated, the uncertainty assigned to each parameter using an alternative 
ht is the absolute difference in its value between the nominal and alternative ht. 

candidates extremely close to the edges of the allowed 


As mentioned in Sect. 4.3 


kinematic region of the Dalitz plot are excluded. The requirement made is that the largest 
of the three |cos(6'ab)| valnes is less than 0.98. A systematic nncertainty dne to this 
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process is estimated by changing the threshold to 0.96, as this excludes a similar additional 
area of the Dalitz plot as the original requirement. 

The systematic uncertainty related to the efficiency model ?^K 7 t) is evaluated 

in four ways. The hrst probes the process by which a smooth curve is produced from 
simulated events; this uncertainty is evaluated using an alternative ht that substitutes the 
non-parametric estimator with a polynomial parameterization. The second uncertainty 
is due to the limited sample size of simulated events. This is evaluated by generating 
several alternative polynomial efficiency models according to the covariance matrix of 
the polynomial model parameters; the spread in parameter values from this ensemble 
is assigned as the uncertainty due to the limited sample size. The third contribution 
is due to possible imperfections in the description of the data by the simulation. This 
uncertainty is assigned using an alternative simultaneous ht that separates the sample into 
three categories according to the year in which the data were collected and the type of Kg 
candidate used. As noted in Sect. the sample recorded during 2011 does not include 
downstream Kg candidates. These sub-samples have different kinematic distributions 
and ^kti) behavior, so this procedure tests the ability of the simulation process 

to reproduce the variation seen in the data. The hnal contribution is due to the re¬ 
weighting procedures used to include the effect of offline selection requirements based on 
information from the RICH detectors, and to correct for discrepancies between data and 
simulation in the reconstruction efficiencies of long and downstream Kg candidates. This 
is evaluated using alternative efficiency models where the relative proportion of the track 
types is altered, and the weights describing the efficiency of selection requirements using 
information from the RICH detectors are modified to account for the limited calibration 
sample size. Additional robustness checks have been performed to probe the description 
of the efficiency function by the simulated events. In these checks the data are divided 
into two equally populated bins of the meson p, px or r] and the amplitude models are 
re-htted using each bin separately. The ht results in each pair of bins are found to be 
compatible within the assigned uncertainties, indicating that the simulated kinematics 
adequately match the data. 

An uncertainty is assigned due to the description of the hardware trigger efficiency in 
simulated events. Because the hardware trigger is not only required to hre on the signal 
decay, it is important that the underlying pp interaction is well described, and a systematic 
uncertainty is assigned due to possible imperfections. This uncertainty is obtained using 
an alternative efficiency model generated from simulated events that have been weighted 
to adjust the fraction where the hardware trigger was hred by the signal candidate. 

The uncertainty due to the description of the combinatorial background is evaluated 
by recomputing the CK 0 K± 7 t=F function using rrijyo sideband events to which 

an alternative kinematic fit has been applied, without a constraint on the D° mass. The 
alternative model is expected to describe the edges of the phase space less accurately, 
while providing an improved description of peaking features. 

An alternative set of models is produced using a threshold of 9 units in the value of 
—2 log C instead of the thresholds of 16 and 25 used for the model building procedure. 
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These models contain more resonances, as fewer are removed dnring the model bnilding 
process. A systematic uncertainty is assigned using these alternative models for those 
parameters which are common between the two sets of models. 

Two parameters of the Flatte dynamical function, which is used to describe the 
ao(980)^ resonance, are hxed to nominal values in the isobar model hts. Alternative 
hts are performed, where these parameters are hxed to different values according to 
their quoted uncertainties, and the largest changes to the ht parameters are assigned as 
systematic uncertainties. 

The effect of resolution in the and coordinates is neglected in the isobar 

model hts, and this is expected to have an effect on the measured K*(892)=*" decay width. 
An uncertainty is calculated using a pseudoexperiment method, and is found to be small. 

The uncertainty due to the yield determination process described in Sect. is measured 
by changing the fractions /m and fc in the isobar model ht according to their statistical 
uncertainties, and taking the largest changes with respect to the nominal result as the 
systematic uncertainty. 

There are two sources of systematic uncertainty due to the amplitude model formalism 
considered. The hrst is that due to varying the meson radius parameters d^o and dehned 
in Sect. 4.f These are changed from d^o = 5.0(GeV/c)“^ 


and dR = 1.5 (GeV/c)“^ to 
2.5 (GeV/c)“^ and 1.0 (GeV/c)“^, respectively. The second is due to the dynamical function 
Tr used to describe the p(1450,1700)^ resonances. These resonances are described by 
the Gounaris-Sakurai functional form in the nominal models, which is replaced with a 
relativistic P-wave Breit-Wigner function to calculate a systematic uncertainty due to this 
choice. 

The uncertainties described above are added in quadrature to produce the total 
systematic uncertainty quoted for the various results. For most quantities the dominant 
systematic uncertainty is due to the meson radius parameters d^o and dR. The largest 
sources of experimental uncertainty relate to the description of the efficiency variation 
across the Dalitz plot. The ht procedure and statistical uncertainty calculation have been 
validated using pseudoexperiments and no bias was found. 

Tables summarizing the various sources of systematic uncertainty and their relative 
contributions are included in Appendix 


5.3 Isobar model results 

The ht results for the best isobar models using the GLASS and LASS parameterizations of the 
Ktt F-wave are given in Tables and Distributions of m‘^^, m'^o^ and are shown 

alongside the best model of the D°—)■ KgK“7r^ mode using the GLASS parameterization in 
Fig-ii In Fig. I^and elsewhere the nomenclature Ri x R 2 denotes interference terms. The 
corresponding distributions showing the best model using the LASS parameterization are 
shown in Fig. Distributions for the D° —)■ KgK+7r“ mode are shown in Figs. and 
Figure shows the GLASS isobar models in two dimensions, and demonstrates that the 
GLASS and LASS choices of Ktt S'-wave parameterization both lead to similar descriptions 
of the overall phase variation. Figures [5]-[^ show distributions distorted by efficiency effects. 
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Figure 5: Distributions of rri^^ (upper left), (upper right) and (lower left) in the 

ROR- 'Ti^ mode with ht curves from the best GLASS model. The solid (blue) curve shows 
the full PDF 4K0K-7t+("^K° 7 t’ while the other curves show the components with the largest 

integrated fractions. 


while Fig. shows the decay rate without distortion. Lookup tables for the complex 
amplitude variation across the Dalitz plot in all four isobar models are available in the 
supplemental material. 

The data are found to favor solutions that have a signihcant neutral Ktt 5-wave 
contribution, even though the exchange (Fig. [^) and penguin annihilation (Fig. [^) 
processes that contribute to the neutral channel are expected to be suppressed. The 
expected suppression is observed for the P-wave K*(892) resonances, with the neutral mode 
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Figure 6: Distributions of (upper left), (upper right) and (lower left) in the 

ROR- 'Ti^ mode with fit curves from the best LASS model. The solid (blue) curve shows the 
full PDF Pkok (m^Ojjj ai^Ti), while the other curves show the components with the largest 
integrated fractions. 


fit fractions substantially lower. The models using the LASS parameterization additionally 
show this pattern for the K*(1410) states. The sums of the ht fractions [^, excluding 
interference terms, in the D° —)■ KgK“7f^ and D° —)■ KgK’''7r“ models are, respectively, 
103% (109%) and 81% (99%) using the GLASS (LASS) Ktt S-wave parameterization. 

Using measurements of the mean strong-phase difference between the D° —)■ KgK^Tt^ 
modes available from i|)(3770) decays 12 , the relative complex amplitudes between each 
resonance in one D° decay mode and its conjugate contribution to the other D° decay 
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Figure 7: Distributions of (upper left), (upper right) and (lower left) in the 

)• KgR"*"?!” mode with ht curves from the best GLASS model. The solid (blue) curve shows 
the full PDF pKgK+Tt- (”^K° 7 t’ ^K 7 t)j while the other curves show the components with the largest 
integrated fractions. 


mode are computed. These values are summarized in Table 

Additional information about the models is listed in Appendices and including 
the interference fractions and decomposition of the systematic uncertainties. The best 
models also contain contributions from the p(1450)^ and p(1700)^ resonances in the 


KgK^ channels, supporting evidence in Ref. 44 of the KK decay modes for these states. 


Alternative models are htted where one p^ contribution is removed from the best models; 
in these the value of —21og£ is found to degrade by at least 162 units. Detailed results 
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Figure 8: Distributions of (upper left), (upper right) and (lower left) in the 

)• KgR"*"?!” mode with fit curves from the best LASS model. The solid (blue) curve shows the 
full PDF pKgK+Tt-("^K^Ti’■''^tiilo tho other curves show the components with the largest 
integrated fractions. 


are tabulated in Appendix 

The Ktt A-wave systems are poorly understood and there is no clear theoretical 
guidance regarding the correct description of these systems in an isobar model. As 
introduced in Sect. 4T, the LASS parameterization is motivated by the Watson theorem, 
but this assumes that three-body interactions are negligible and is not, therefore, expected 
to be precisely obeyed in nature. The isobar models using the GLASS parameterization favor 
solutions with qualitatively similar phase behavior to those using the LASS parameterization. 
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Figure 9: Decay rate and phase variation across the Dalitz plot. The top row shows 
|A^K0K±7tT ^K 7 i)P ™ tiest GLASS isobar models, the center row shows the phase 

behavior of the same models and the bottom row shows the same function subtracted from the 
phase behavior in the best LASS isobar models. The left column shows the D^—)• KgK~7t^ mode 
with )• KgK'''7r“ on the right. The small inhomogeneities that are visible in the bottom row 
relate to the GLASS and LASS models preferring slightly different values of the K*(892)^ mass 
and width. 
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Table 7: Modulus and phase of the relative amplitudes between resonances that appear in both 


the K°K-7t+ and K^K+Tf 

of measured in 4)(3770) decays 

statistical uncertainty. The first uncertainties are statistical and the second systematic. 


modes. Relative phases are calculated using the value 
and the uncertainty on this value is included in the 


12 


Relative 


amplitude 


GLASS 

LASS 

^(K*(892)-) 

.4(K*(892)+) 

mod 
arg (°) 

0.582 ±0.007 ±0.008 
-2 ±15 ±2 

0.576 ±0.005 ±0.010 
-2±15 ± 1 

.4(K*(1410)-) 
^{K*(1410) + ) 

mod 
arg (°) 

0.64 ±0.08 ±0.22 

52 ±17 ±20 

0.90 ±0.08 ±0.15 

62 ± 16 ± 6 

-4((Kg")s-wave) 

-4((Kg")J_wave) 

mod 
arg (°) 

0.54 ±0.06 ±0.26 
-100 ±20 ±40 

0.59 ±0.05 ±0.08 
-44 ± 17 ± 10 

.4(K* (892)0) 

^(K* (892)0) 

mod 
arg (°) 

1.12±0.05 ±0.11 

-78 ±16 ±10 

1.17±0.04±0.05 
-75 ± 15 ± 2 

.4(K* (1410)0) 

.4(K* (1410)0) 

mod 
arg (°) 

0.60 ±0.05 ±0.12 
-9±16 ± 14 

0.62 ±0.09 ±0.12 
-23 ± 17 ± 11 

^(K* (1430)0) 

^(k; (1430)0) 

mod 
arg (°) 

1.1±0.1±0.5 

31 ± 17 ± 12 

— 

^((K+7t“)s_„ave) 

^({K 7t”*”)5_'wave) 

mod 
arg (°) 

0.87±0.08±0.14 

49 ± 25 ± 16 

0.78 ±0.06 ±0.18 

68 ± 16 ± 6 

^(ao(980) + ) 
.4(ao(980)-) 

mod 
arg (°) 

— 

2.1±0.2±0.6 

42 ± 16 ± 5 

.4(ao(1450)+) 

^(ao(1450)-) 

mod 
arg (°) 

0.49 ±0.06 ±0.28 
-60 ±19 ±34 

1.14 ±0.16 ±0.30 

-63 ±20 ± 19 

^(p(1450) + ) 
.4(p(1450)-) 

mod 
arg (°) 

0.86 ±0.16 ±0.26 

110 ±20 ±50 

— 

^(p(1700) + ) 
.4(p(1700)-) 

mod 
arg (°) 

1.6 ±0.4 ±0.4 

70 ± 20 ± 70 

— 


This is illustrated in Fig. 10, which also shows the GLASS forms obtained in hts to 


D°—)■ KgTf^Tt decays by the BaBar collaboration 37138 and previously used in fits to the 


D° —)■ KgK^Tt^ decay modes 12 . This hgure shows that the GLASS functional form has 


substantial freedom to produce different phase behavior to the LASS form, but that this is 
not strongly favored in the D° —)■ KgK^Tt^ decays. The good quality of fit obtained using 
the LASS parameterization indicates that large differences in phase behavior with respect 
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Figure 10: Comparison of the phase behavior of the various Ktt 5-wave parameterizations used. 
The solid (red) curve shows the LASS parameterization, while the dashed (blue) and dash-dotted 
(green) curves show, respectively, the GLASS functional form fitted to the charged and neutral 
5-wave channels. The final two curves show the GLASS forms fitted to the charged Ktt 5-wave 
in —)■ KgTT+TT” decays in Ref. 37 (triangular markers, purple) and Ref. 38 (dotted curve, 
black). The latter of these was used in the analysis of —)• KgK^Tt^ decays by the CLEO 
collaboration 


12 


to Ktt scattering data 45,46 are not required in order to describe the D° —)■ KgK^Tt^ 
decays. A similar conclusion was drawn in Ref. for the decay D’*'—)• K“7r’''7r+, while 
Ref. 58 found behavior inconsistent with scattering data using the same decay mode 
but a slightly different technique. Ref. 59 studied the Ktt 5-wave in x~ —)■ KgTfv,- decays 


and found that a parameterization based on the LASS Ktt scattering data, but without a 
real production form factor, gave a poor description of the x~ decay data. 

The quality of ht for each model is quantihed by calculating using a dynamic 
binning scheme. The values are summarized in Table while the binning scheme and 
two-dimensional quality of ht are shown in Appendix]^ This binning scheme is generated 
by iteratively sub-dividing the Dalitz plot to produce new bins of approximately equal 
population until further sub-division would result in a bin population of fewer than 15 
candidates, or a bin dimension smaller than 0.02 GeVyc^ in or This minimum 

size corresponds to hve times the average resolution in these variables. 

The overall ht quality is slightly better in the isobar models using the GLASS Ktt 5-wave 
parameterization, but this is not a signihcant effect and it should be noted that these 
models contain more degrees of freedom, with 57 parameters htted in the hnal GLASS 
model compared to 50 when using the LASS parameterization. 
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Table 8: Values of x^/bin indicating the fit quality obtained using both Krt S'-wave parameteri- 
zations in the two decay modes. The binning scheme for the )> KgK^TT"^ (D®—>■ KgK'^Tt”) 
mode contains 2191 (2573) bins. 


Isobar model 
GLASS LASS 

D0^K0K-7r+ 1.12 1.10 
DO^KOK+tt- 1.07 1.09 

6 Additional measurements 

In this section, several additional resnlts, inclnding those derived from the amplitnde 
models, are presented. 

6.1 Ratio of branching fractions measurement 

The ratio of branching fractions 


_ KOK+tT 

= 


K0K-7r+)’ 


( 18 ) 


and the restricted region ratio defined in the same region near the K*(892)^ resonance 

as the coherence factor i7K*K (Sect. 4.2), are also measured. The efficiency correction due 

is evaluated using the best isobar models. 


to the reconstruction efficiency 
and the difference between the results obtained with the two Ktt S'-wave parameterizations 
is taken as a systematic uncertainty in addition to those effects described in Sect. |5.2[ This 
efficiency correction modifies the ratio of yields quoted in Table [T] by approximately 3%. 

The two ratios are measured to be 


'^KgKTt = 0.655 ± 0.004 (stat) ± 0.006 (syst), 
'Bk*k = 0.370 ± 0.003 (stat) ± 0.012 (syst). 

These are the most precise measurements to date. 


6.2 Coherence factor and CP-even fraction results 


The amplitude models are used to calculate the coherence factors i^KgKn Rk*Ki and the 


strong-phase difference ^KgKTi ~ <^k*k 


as described in Sect. 


4.2 


The results are summarized 


in Table alongside the corresponding values measured in if) (3770) decays by the CLEO 
collaboration. Lower, but compatible, coherence is calculated using the isobar models than 
was measured at CLEO, with the discrepancy larger for the coherence factor calculated over 
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Table 9: Coherence factor observables to which the isobar models are sensitive. The third column 
summarizes the CLEO results measured in quantum-correlated decays [^, where the uncertainty 
^KgKTi “ '^K*K is calculated assuming maximal correlation between (ixgKTt 5 k*k- 


Variable GLASS LASS CLEO 


^KgKTt 0.573 ±0.007 ±0.019 

Rk*k 0.831 ±0.004 ±0.010 

5KOKn-<5K*K (0.2 ±0.6 ±1.1)° 


0.571 ±0.005 ±0.019 
0.835 ±0.003 ±0.011 
(-0.0 ±0.5 ±0.7)° 


0.73±0.08 

1.00±0.16 

(-18±31)° 


the full phase space. The results from the GLASS and LASS isobar models are very similar, 
showing that the coherence variables are not sensitive to the Ktt S'-wave parameterization. 

The coherence factor the ratio of branching fractions B^oku^ combined 

with the mean phase difference between the two hnal states measured in i|)(3770) decays 
to calculate the CP-even fraction P+, dehned in Eq. 16, which is determined to be 


12 


P+ = 0.777 ± 0.003 (stat) ± 0.009 (syst), 

using the GLASS amplitude models. A consistent result is obtained using the alternative 
(lass) amplitude models. This model-dependent value is compatible with the direct 
measurement using only i|)(3770) decay data 


12,51 


6.3 SU(3) flavor symmetry tests 

SU(3) flavor symmetry can be used to relate decay amplitudes in several D meson decays, 
such that a global ht to many such amplitudes can provide predictions for the neutral 
and charged K*(892) complex amplitudes in D°—;• KgK^Tt"^ decays |^[^. Predictions are 
available for the K*+K“, K*“K+, and complex amplitudes, where K* refers 

to the K*(892) resonances. There are therefore three relative amplitudes and two relative 
phases that can be determined from the isobar models, with an additional relative phase 
accessible if the isobar results are combined with the CLEO measurement of the mean 
strong phase difference 


12 


The results are summarized in Table 10 


The isobar model results are found to follow broadly the patterns predicted by SU(3) 
flavor symmetry. The amplitude ratio between the K*(892)’'' and K (892)° resonances, 
which is derived from the D°—)■ KgK^Tr’*' isobar model alone, shows good agreement. The 
two other amplitude ratios additionally depend on the ratio i^KgRTu ^ind these are more 
discrepant with the SU(3) predictions. The relative phase between the charged and neutral 
K* (892) resonances shows better agreement with the flavor symmetry prediction in the 
D°—> KgK+7r“ mode, where both resonances have clear peaks in the data. The GLASS and 
LASS isobar models are found to agree well, suggesting the problems are not related to the 
Ktt S'-wave. 
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Table 10; SU(3) flavor symmetry predictions and results. The uncertainties on phase difference 
predictions are calculated from the quoted magnitude and phase uncertainties. Note that some 
theoretical predictions depend on the ri-r|' mixing angle 6*T,-n' and are quoted for two different 
values. The bottom entry in the table relies on the CLEO measurement |12| of the coherence 
factor phase and the uncertainty on this phase is included in the statistical uncertainty, 

while the other entries are calculated directly from the isobar models and relative branching 
ratio. Where two uncertainties are quoted the first is statistical and the second systematic. 


Theory Experiment 


Ratio 

9^.^, = 19.5° 

9^.^, = 11.7° 

GLASS 

LASS 

|^(K*(892)-K+) 

|.4(K*(892)+K-) 

0.685 ±0.032 

0.685 ±0.032 

0.582 ±0.007 ±0.007 

0.576 ±0.005 ±0.010 

|yi(K*(892)0K0) 
|y^l(K*(892) + K-)| 

0.138 ±0.033 

0.307 ±0.035 

0.297 ±0.010 ±0.024 

0.295 ±0.009 ±0.014 

|y!l(K*(892)0K0) 

|.4(K*(892)+K-)t 

0.138 ±0.033 

0.307 ±0.035 

0.333 ±0.008 ±0.016 

0.345 ±0.007 ±0.010 

Argument 

Theory (°) 

Experiment (°) 

.4(K*(892)0k0) 

.4(K*(892)+K-) 

151 ±14 

112±8 

72 ± 2 ± 4 

78.5 ± 2.0 ± 2.8 

.4(K*(892)0k0) 

.4(K*(892)-K+) 

-9 ±13 

-37 ±6 

-4±2±9 

5.0±1.7± 1.4 

.4(K*(892)0k0) 

.4(K*(892)0k0) 

180 

180 

-78 ± 16 ± 10 

-75 ±15 ±2 


6.4 CP violation tests 


Searches for time-integrated CP-violating effects in the resonant structure of these decays 
are performed using the best isobar models. The resonance amplitude and phase parameters 
Or and 0R are substituted with aR(l ± Aor) and 0 r ± A0r, respectively, where the signs 
are set by the flavor tag. The convention adopted is that a positive sign produces the D° 
complex amplitude. The full £t results are tabulated in Appendix [Pj 

A subset of the A parameters is used to perform a test against the no-CP violation 
hypothesis: only those parameters corresponding to resonances that are present in the best 
isobar models using both the GLASS and LASS Ktt S-wave parameterizations are included. 
The absolute difference |Aglass ~ AlassI is assigned as the systematic uncertainty dne 
to dependence on the choice of isobar model. This snbset of parameters is shown in 
Table IT where the change in ht fraction between the D° and solutions is included for 
illustrative purposes. In the test the statistical and systematic uncertainties are added 
in qnadrature. 

Using the best GLASS (LASS) isobar models the test result is y^/ndf = 30.5/32 = 0.95 
(32.3/32 = 1.01), corresponding to a p-valne of 0.54 (0.45). Therefore, the data are 
compatible with the hypothesis of CP-conservation. 
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Table 11: CP violation fit results. Results are only shown for those resonances that appear in both the GLASS and LASS models. The 
first uncertainties are statistical and the second systematic; the only systematic uncertainty is that due to the choice of isobar model. 
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ao(980)+ 0.06±0.04 ± 0.01 0.052±0.025 ± 0.008 -3±5±2 -0.9±3.1±2.2 2.2±2.8 ±2.4 4.6±3.3 ±2.4 

ao(1450)+ -0.11 ±0.10 ± 0.04 -0.07± 0.07 ± 0.04 10±8±5 5±6±5 -0.21 ± 0.30 ± 0.23 -0.4±0.4±0.2 

p(1700)+ -0.03± 0.13 ±0.09 -0.12± 0.13 ± 0.09 4±6±2 2±5±2 -0.07±0.25± 0.19 -0.27± 0.27 ± 0.19 



7 Conclusions 


The decay modes D° —)■ KgK^Tt^ have been studied using unbinned, time-integrated, hts 
to a high purity sample of 189 670 candidates, and two amplitude models have been 
constructed for each decay mode. These models are compared to data in a large number 
of bins in the relevant Dalitz plots and a fest indicates a good description of the data. 

Models are presented using two different parameterizations of the Ktt S'-wave systems, 
which have been found to be an important component of these decays. These systems are 
poorly understood, and comparisons have been made to previous results and alternative 
parameterizations, but the treatment of the Ktt S'-wave is found to have little impact on 
the other results presented in this paper. The large fractions attributed to the neutral 
Ktt S'-wave channels could indicate larger than expected contributions from the penguin 
annihilation diagrams shown in Fig. [^. 

The models are seen to favor small, but signihcant, contributions from the 
p(1450,1700)^ —)■ KgK^ resonances, modes which were seen by the OBELIX experi¬ 
ment but are not well established. All models contain clear contributions from 
both the K*(892)=*" and K*(892)° resonances, with the K*(892)° contribution found to be 
suppressed as expected from the diagrams shown in Fig. This allows the full set of 
amplitudes in these decays that are predicted by SU(3) flavor symmetry to be tested, in 
contrast to the previous analysis by the CLEO collaboration [^. Partial agreement is 
found with these predictions. 

The ratio of branching fractions between the two D° —)■ KgK^TrfF modes is also measured, 
both across the full Dalitz plot area and in a restricted region near the K*(892)^ resonance, 
with much improved precision compared to previous results. 

Values for the D° —)■ KgKyt coherence factor are computed using the amplitude models, 
again both for the whole Dalitz plot area and in the restricted region, and are found to be 
in reasonable agreement with direct measurements by CLEO 12 using quantum-correlated 


r[)(3770) —)■ D°D° decays. The CP-even fraction of the D°—)■ KgKTr decays is also computed, 
using input from the quantum-correlated decays, and is found to be in agreement with the 
direct measurement 12,51 • A search for time-integrated CP violation is carried out using 


the amplitude models, but no evidence is found with either choice of parameterization for 
the Ktt S'-wave. 

The models presented here will be useful for future D°-D° mixing, indirect CP violation 
and CKM angle 7 studies, where knowledge of the strong-phase variation across the Dalitz 
plot can improve the attainable precision. These improvements will be particularly valuable 
for studies of the large dataset that is expected to be accumulated in Run 2 of the LHC. 


Acknowledgements 

We express our gratitude to our colleagues in the CERN accelerator departments for the 
excellent performance of the LHC. We thank the technical and administrative staff at the 
LHCb institutes. We acknowledge support from CERN and from the national agencies: 
CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 (France); 


29 






BMBF, DFG and MPG (Germany); INFN (Italy); FOM and NWO (The Netherlands); 
MNiSW and NGN (Poland); MEN/IFA (Romania); MinES and FANG (Russia); MinECo 
(Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (United Kingdom); NSF 
(USA). We acknowledge the computing resources that are provided by GERN, IN2P3 
(France), KIT and DESY (Germany), INFN (Italy), SURF (The Netherlands), PIC (Spain), 
GridPP (United Kingdom), RRCKI (Russia), CSCS (Switzerland), IFIN-HH (Romania), 
CBPF (Brazil), PL-GRID (Poland) and OSC (USA). We are indebted to the communities 
behind the multiple open source software packages on which we depend. We are also 
thankful for the computing resources and the access to software R&D tools provided by 
Yandex LEG (Russia). Individual groups or members have received support from AvH 
Foundation (Germany), EPLANET, Marie Sklodowska-Curie Actions and ERG (European 
Union), Conseil General de Haute-Savoie, Labex ENIGMASS and OCEVU, Region 
Auvergne (France), RFBR (Russia), XuntaGal and GENCAT (Spain), The Royal Society 
and Royal Commission for the Exhibition of 1851 (United Kingdom). We acknowledge 
the use of the Advanced Research Computing (ARC) facilities (Oxford) in carrying out 
this work. 


30 



Appendices 


A Additional isobar formalism information 


This appendix contains Table 12, which summarizes the nominal values used for the various 
resonance and form factor parameters. 


B Additional isobar model information 


This appendix contains additional information about the various isobar model parameters 
that are used and allowed to vary freely in the model fits, e.g. resonance mass and width 
parameter values, and parameters of the GLASS and LASS Ktt S'-wave functional forms. 
Tables 13 and 14 summarize the most significant interference terms in the D° —)• 


and D° —)• KgK+Tt models, respectively. Table 15 defines the matrices U used to define 


the LASS Ktt S'-wave form factor. Tables 16 (GLASS) and 17 (LASS) summarize the various 
resonance and form factor parameters. The nominal values that are used in Gaussian 
constraint terms are given in Appendix [A| 


Figure pT] shows the smooth functions that describe the combinatorial background in 
the isobar model fits. Figure illustrates the two-dimensional quality of fit achieved in 
the four isobar models and shows the binning scheme used to derive y^/bin values. 

The changes in —21og£ obtained in alternative models where one p contribution is 


removed are given in Table 18 
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Figure 11: Smooth functions, used to describe the combinatorial back¬ 

ground component in the —?■ KgK^Tf’' (left) and —)> KgK+Tf (right) amplitude model 
fits. 
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Table 12: Nominal values for isobar model parameters that are fixed in the model fits, or used 
constraint terms. These values are taken from Refs. [3|[35|[44|[4^ as described in Sect.|^ 


Parameter Value 


K*(892)± 

Tr 

891.66±0.26 

50.8±0.9 

MeV/c2 

MeV/c2 

K*(1410)± 

"iR 

Tr 

1.414±0.015 

0.232±0.021 

GeV/c2 

GeV/c2 


^R 

Tr 

1.435±0.005 

0.279±0.006 

GeV/c2 

GeV/c2 

K* (892)0 

"iR 

Tr 

895.94±0.22 

48.7±0.8 

MeV/c2 

MeV/c2 

K*(1410)0 

"iR 

Tr 

1.414±0.015 

0.232±0.021 

GeV/c2 

GeV/c^ 

K*(1430)0 

^R 

Tr 

1.4324±0.0013 

0.109±0.005 

GeV/c2 

GeV/c^ 


"IR 

Fr 

1.435±0.005 

0.279±0.006 

GeV/c^ 

GeV/c2 

Ktt S'-wave 

r 

a 

1.8±0.4 

1.95±0.09 

(GeV/c)-^ 

(GeV/c)-i 

ao(980)± 

m-R 

9r]Ti 

0.980±0.020 

324±15 

GeV/c2 

MeV 


^KK 

9r\7i 

1.03±0.14 


a2(1320)± 

m-R 

Fr 

1.3181±0.0007 

0.1098±0.0024 

GeV/c^ 

GeV/c2 

ao(1450)± 

rriR 

Fr 

1.474±0.019 

0.265±0.013 

GeV/c^ 

GeV/c2 

p(1450)± 

rriR 

Fr 

1.182±0.030 

0.389±0.020 

GeV/c2 

GeV/c2 

p(1700)± 

rriR 

Fr 

1.594±0.020 

0.259±0.020 

GeV/c2 

GeV/c2 


in 
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Table 13: Interference fractions for the —?■ KgK Tf*" mode. The first uncertainties are statistical and the second systematic. Only 

the 25 largest terms are shown. 
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Table 14: Interference fractions for the —?■ KgR"*"?! mode. The first uncertainties are statistical and the second systematic. Only 

the 25 largest terms are shown. 
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Table 15: Matrices U relating the fit coordinates b' to the LASS form factor coordinates b = Ub' 
defined in Sect. 


/-0.460 0.702 -0.543\ /-0.452 -0.676 -0.582\ 

0.776 0.197 -0.631 0.776 0.0243 -0.631 

\-0.433 -0.711 -0.554/ \-0.440 0.737 -0.513/ 




Figure 12: Two-dimensional quality-of-fit distributions illustrating the dynamic binning scheme 
used to evaluate y^. The variable shown is where di and pi are the number of events and 

yPi 

the fitted value, respectively, in bin i. The (D^—)• KgR"*"?!”) mode is shown in 

the left (right) column, and the GLASS (LASS) isobar models are shown in the top (bottom) row. 
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Table 16: Additional fit parameters for GLASS models. This table does not include parameters 
that are fixed to their nominal values. The hrst uncertainties are statistical and the second 
systematic. 


Parameter Value 


K*(892)± 

"iR 

Tr 

893.1 ±0.1 ±0.9 

46.9 ±0.3 ±2.5 

MeV/c2 

MeV/c2 

K*(1410)± 

Tr 

210 ±20 ±60 

MeV/c2 


F 

1.785 (fixed) 



a 

4.7±0.4±1.0 

(GeV/c)-^ 

(Ks7T)|„ave 

4>f 

0.28 ±0.05 ±0.19 

rad 


4>s 

2.8 ±0.2 ±0.5 

rad 


r 

-5.3 ±0.4 ±1.9 

(GeV/c)-^ 

K*(1410)° 

Tr 

1426 ± 8 ± 24 

270 ± 20 ± 40 

MeV/c2 

MeV/c2 


F 

0.15 ±0.03 ±0.14 



a 

4.2 ±0.3 ±2.8 

(GeV/c)-^ 

(K7r)°__ 

4>f 

-2.5±0.2 ± 1.0 

rad 


4>s 

-1.1±0.6 ± 1.3 

rad 


r 

-3.0±0.4±1.7 

(GeV/c)-^ 

ao(1450)± 


1430 ±10 ±40 

MeV/c2 

p(1450)± 

Tr 

410 ±19 ±35 

MeV/c2 

p(1700)± 

»^R 

1530 ±10 ±40 

MeV/c2 
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Table 17: Additional fit parameters for LASS models. This table does not include parameters 
that are fixed to their nominal values. The first uncertainties are statistical and the second 
systematic. 


Parameter Value 


K*(892)± 

Tr 

893.4 ±0.1 ± 1.1 
47.4 ±0.3 ±2.0 

MeV/c2 

MeV/c2 

K*(1410)^ 

mji 

1437±8± 16 

MeV/c2 


b[ 

60 ± 30 ± 40 



b'2 

4±1±5 




3.0±0.2±0.7 


K*(1410)° 

mji 

1404 ± 9 ± 22 

MeV/c2 


K 

130 ±30 ±80 



b'2 

-6±1 ± 14 



&3 

2.5±0.1 ± 1.4 


Ktt 5'-wave 

r 

1.2±0.3±0.4 

(GeV/c)-i 

ao(980)± 

rriK 

925±5±8 

MeV/c2 

ao(1450)± 

rriK 

Tr 

1458 ± 14 ± 15 

282 ±12 ± 13 

MeV/c2 

MeV/c2 

p(1450)± 

^R 

1208±8±9 

MeV/c2 

p(1700)± 

^R 

1552 ±13 ±26 

MeV/c2 


Table 18: Change in —21og£ value when removing a p resonance from one of the models. 


Ktt S'-wave Removed 

parameterization Decay mode resonance A(—21og£) 


LASS 


GLASS 


D°^K0K-7t+ p(1450)- 338 

D°^K0K+7t- p(1700)+ 235 


K0K-7t+ 

K0K-7t+ 


p(1450)- 238 

p(1700)- 162 

p(1450)+ 175 

p(1700)+ 233 
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C Systematic uncertainty tables 

This appendix includes tables summarizing the various contributions to the systematic 
uncertainties assigned to the various results. The table headings correspond to the 
uncertainties discussed in Sect. 15.21 with some abbreviations to allow the tables to be 
typeset compactly. Dehnitions of the various abbreviations are given in Table The 


quantity ‘DFF’ listed in the tables is the sum of fit fractions from the various resonances, 
excluding interference terms. Tables 20 (GLASS) and (LASS) show the results for the 
complex amplitudes and ht fractions in the D° —)• models, Tables 22 (GLASS) 

(lass) show the corresponding values for the D°—)• KgK+7r“ models and Tables 24 


23 


25 summarize the uncertainties for the parameters that are not specihc to a decay 


and 
and 
mode. 

In each of these tables the parameter in question is listed on the left, followed by the 
central value and the corresponding statistical (first) and systematic (second) uncertainty. 
The subsequent columns list the contributions to this systematic uncertainty, and are 
approximately ordered in decreasing order of significance from left to right. 


D CP violation fit results 


This appendix contains Table 26 which summarizes the full fit results of the CP violation 
searches described in Sect. 16.41 
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Table 19; Listing of abbreviations required to typeset the systematic uncertainty tables. 


Abbreviation 
max(| cos |) 

Efficiency 

Joint 

Weights 

Comb. 

—2 log£ 

Flatte 

/m,/c 

(Jdo, (Jr 


Description 

Variation of the cut that excludes the boundary regions of the Dalitz 
plot. 

Two efficiency modelling uncertainties added in quadrature; using an 
alternative parameterization, and accounting for the limited size of the 
simulated event sample. 

Uncertainty obtained by simultaneously htting disjoint sub-sets of the 
dataset, separated by the year of data-taking and type of Kg daughter 
track, with distinct efficiency models. 

Three uncertainties related to the re-weighting of simulated events used 
to generate the efficiency model added in quadrature. 

These account for; incorrect simulation of the underlying pp interaction, 
uncertainty in the relative yield of long and downstream Kg candidates, 
and uncertainty in the efficiency of selection requirements using informa¬ 
tion from the RICH detectors. 

Using an alternative combinatorial background model. 

Using a more complex alternative model where the threshold in 
A(—21og£) for a resonance to be retained is reduced to 9 units. 
Variation of the Flatte lineshape parameters for the ao(980)^ resonance 
according to their nominal uncertainties. 

Variation of the mistag and combinatorial background rates according 
to their uncertainties in the mass ht. 

Variation of the meson radius parameters. 

Switching to a Breit-Wigner dynamical function to describe the 
p(1450,1700)^ resonances. 
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Table 20: Systematic uncertainties for complex amplitudes and fit fractions in the —)• KgK 7T+ model using the GLASS parameteri¬ 

zation. The headings are defined in Table PI 
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Table 21: Systematic uncertainties for complex amplitudes and fit fractions in the )• KgK 7T+ model using the LASS parameteri¬ 
zation. The headings are defined in Table PI 
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Table 22: Systematic uncertainties for complex amplitudes and fit fractions in the —)• KgR"*"?! model using the GLASS parameteri¬ 

zation. The headings are defined in Table PI 
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OR 1.04 ±0.12 ±0.32 0.23 0.09 0.07 0.15 0.09 0.05 0.06 0.02 0.01 0.01 

p(1700)+ 4 > k [°) 4 ±11 ±20 4.6 13.9 7.3 2.4 8.4 6.9 2.2 1.6 1.8 2.5 

FF [%] 0.5±0.1±0.5 0.39 0.12 0.07 0.13 0.02 0.06 0.05 0.02 0.01 0.11 

xVbin 1-07 1.09 1.07 1.06 1.04 1.06 - - 1.07 1.07 1.12 

DFF [%] 80.7 78.5 71.3 84.9 83.4 82.9 - - 79.9 80.5 81.5 












































Table 23: Systematic uncertainties for complex amplitudes and fit fractions in the > KgK+Tt model using the LASS parameteri¬ 
zation. The headings are defined in Table PI 
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OR 1.25T0.15T 0.33 0.22 0.12 0.07 0.12 0.12 0.02 0.02 0.04 0.11 0.02 

p(1700)+ 0R (°) 39±9± 15 9.5 8.4 5.1 4.0 2.0 1.3 0.1 0.4 2.1 1.1 

FF [%] 0.53T0.11 ±0.23 0.15 0.08 0.06 0.07 0.10 0.01 0.07 0.03 0.02 0.01 

xVbin 1.09 1.10 1.07 1.08 - - 1.09 1.14 1.09 1.09 1.09 

DFF [%] 99.0 99.1 104.8 95.6 - - 98.9 100.7 99.3 99.3 99.4 


















Table 24: Systematic uncertainties for shared parameters, coherence and relative branching ratio observables in the GLASS models. 
The systematic uncertainty on the K*(892)^ width due to neglecting resolution effects in the nominal models is 0.5MeV/c^. 
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Table 25: Systematic uncertainties for shared parameters, coherence and relative branching ratio observables in the LASS models. 
The systematic uncertainty on the K*(892)^ width due to neglecting resolution effects in the nominal models is 0.6MeV/c^. 
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Table 26: Full CP violation lit results as described in Sect. 6.4 The only uncertainties included are statistical. 
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Supplemental material 

This is divided into two parts: lookup tables for the complex amplitude and covariance 
information, each for the four quoted amplitude models. These are available at Ref. j^. 
No correlation information is included for systematic uncertainties. 


Isobar model lookup tables 


The lookup table filenames are listed in Table ^ 
file glass_fav_lookup.txt are: 


As an example, the first five lines of the 


# S-wave: GLASS, mode: D0->KSK-pi+ (FAV) 

# mDO = 1.86486; mKS = 0.497614; mK = 0.493677; mPi = 0.13957018 GeV/c''2 

# m~2(Kpi) GeV~2/c~4, m~2(KSpi) GeV"2/c~4, |amp|''2 arb. units, arg(amp) rad 
0.300625,0.300625,0.OOOOOOe+00,0.000000 

0.300625,0.301875,0.OOOOOOe+00,0.000000 

The first three lines are comments, describing which D° decay mode and isobar model this 
file corresponds to, giving the precise nominal masses used in the fit and, finally, defining 
the data fields in the remainder of the file. As this shows, the models are evaluated on a 
grid with a spacing of 0.00125 GeV^/c^. 


Covariance information 

A reduced covariance matrix is presented for each isobar model, tabulating the correlations 
between the complex amplitudes These are listed in hies named analogously to 

those in Table e.g. glass_fav_covariance.txt. An example hrst four lines: 

# S-wave: GLASS, mode: D0->KSK-pi+ (FAV) 

# X , y , cov(x,y) 

K(0)*(1430)+_Amp,K(0)*(1430)+_Amp,2.195e-03 
K(0)*(1430)+_Amp,K(0)*(1430)+_Phase,1.147e-01 

i.e. a similar format to the lookup tables. Note that the Ktt A-wave contributions are 
tabulated as K(0)*(1430)+ and K(0) * (1430)bar0. 


Table 27: Lookup table filenames. 


Ktt A-wave parameterization 
D° decay mode GLASS LASS 

D° —)■ KgK^Tf*" glass_fav_lookup.txt lass_fav_lookup.txt 
D° —)■ KgK+Tt” glass_sup_lookup.txt lass_sup_lookup.txt 
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